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A DIMENSIONAL STUDY OF DISK GALAXIES 
(AN EXCUSE TO TALK ABOUT THE HUBBLE SEQUENCE) 

X. Hernandez,^ B. Cervantes-Sodi ^ 
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RESUMEN 

El resumen sera traducido al espanol por los editores. We present a highly simphfied model of the 
dynamical structure of a disk galaxy where only two parmeters fully determine the solution, mass and angular 
momentum. We show through simple physical scalings that once the mass has been fixed, the angular momen- 
tum parameter A is expected to regulate such critical galactic disk properties as colour, thickness of the disk 
and disk to bulge ratio. It is hence expected to be the determinant physical ingeredient resulting in a given 
Hubble type. A simple analytic estimate of A for an observed system is provided. An explicit comparison of the 
distribution of several galactic parameters against both Hubble type and A is performed using observed galax- 
ies. Both such distributions exhibit highly similar characteristics for all galactic properties studied, suggesting 
A as a physically motivated classification parameter for disk galaxies. 

ABSTRACT 

We present a highly simplified model of the dynamical structure of a disk galaxy where only two parmeters fully 
determine the solution, mass and angular momentum. We show through simple physical scalings that once 
the mass has been fixed, the angular momentum parameter A is expected to regulate such critical galactic disk 
properties as colour, thickness of the disk and disk to bulge ratio. It is hence expected to be the determinant 
physical ingeredient resulting in a given Hubble type. A simple analytic estimate of A for an observed system 
is provided. An explicit comparison of the distribution of several galactic parameters against both Hubble type 
and A is performed using observed galaxies. Both such distributions exhibit highly similar characteristics for 
all galactic properties studied, suggesting A as a physically motivated classification parameter for disk galaxies. 
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1. INTRODUCTION 

The first thing an astronomer wants to know 
about a galaxy is its morphological type, typically 
expressed through the Hubble classification scheme, 
introduced about 80 years ago. 

The continued usage of the Hubble classification 
owes its success to the fact that on hearing the Hub- 
ble entry of a galaxy, one immediately forms an im- 
age of the type of galaxy that is being talked about. 
A wide variety of physical features of galaxies show 
good monotonic correlations with Hubble type, de- 
spite the presence of significant overlap and disper- 
sion. To mention only a few, total magnitudes de- 
crease towards later types and colours become bluer 
while gas fractions diminish from early to late types 
(e.g. Roberts & Haynes 1994). Bulge magnitudes 
and bulge to disk ratios decrease (e.g. Pahre et al. 

2004) and disks become thinner (e.g. Kregel et al. 

2005) in going to later types. 

However, galactic classification schemes suffer 
from several shortcomings. First, the diagnostics 
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which enter into the assignment of type, do so in a 
fundamentally subjective manner. One has to take 
the images to an expert who will somehow integrate 
various aspects of the galaxy in question to arrive 
at the type parameter. That different authors gen- 
erally agree, together with the correlations between 
type and various galactic parameters is indicative 
of a solid physical substrate to galaxy classification 
schemes, the particular nature of which however, has 
remained elusive. Second, type parameters are es- 
sentially qualitative which make it difficult to relate 
type to definitive quantitative aspects of a galac- 
tic system, and bring into question the validity of 
any statistical mathematical analysis performed on 
galactic populations based on type. Finally, the type 
given to a galaxy is highly dependent on the infor- 
mation one has regarding the system. This stems 
from the inputs that determine the galactic type, 
the relevance of the bulge, typically redder than the 
disk, means that when observing at longer rest frame 
bands galaxies shift to earlier types, an effect com- 
pounded with the high wavelength sensitivity of the 
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bright HII regions which define the structure and 
morphology of the spiral arms. Similarly, projec- 
tion effects become important, as the structure of the 
arms completely disappears in edge-on disks, elimi- 
nating one of the diagnostics relevant to the assign- 
ment of type. 

Monotonic trends with Hubble type suggest that 
after having fixed the mass, there might exist one 
other parameter whose variation gives rise to the 
Hubble sequence. Many approaches to galaxy for- 
mation and evolution have appeared over the years, 
generally reaching the conclusion that it is the an- 
gular momentum of a galactic system what deter- 
mines its main characteristics, e.g. Fall & Efstathiou 
(1980), Firmani et al. (1996). 

2. THEORETICAL FRAMEWORK 

Having chosen the angular momentum as our sec- 
ond parameter wc now construct a simple physical 
model for a galaxy which allows to estimate this pa- 
rameter directly from observed galactic properties. 

In the interest of capturing the most essential 
physics of the problem, the model we shall develop 
will be the simplest one could possibly propose, little 
more than a dimensional analysis of the problem. 
We shall model only two galactic components, the 
first one a disk having a surface mass density E(r) 
satisfying: 



(1) 



Where r is a radial coordinate and Eq a-nd Rd are 
two constants which are allowed to vary from galaxy 
to galaxy. The second component will be a fixed 
dark matter halo having an isothermal p{r) density 
profile, and responsible for establishing a rigorously 
flat rotation curve Vd throughout the entire galaxy, 
such that 
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and a halo mass profile M{r) oc r. Wc shall refer 
not to the total angular momentum L, but to the 
dimensionless angular momentum parameter 



A = 



L I E |i/2 



(3) 



where E is the total energy of the configuration. 
We now assume | i? | is given by 1/2 the potential 
energy of the halo, which dominates through a halo 
to disk mass ratio Mn/Md = ~ 25. If the specific 
angular momenta of disk and halo are equal. Id = 
Ih, as would be the case for an initially well mixed 



protogalactic state, we can introduce a disk TuUy- 
Fisher relation: Md = ATFV^^^ into to yield: 



A = 21. 



Rd/kpc 



[Vd/kms-^fl^' 



(4) 



Where Atf and F were calibrated from Galac- 
tic parameters. The above equation allows a direct 
estimate of A, a dimensionless numerical parameter 
with a clear physical interpretation for any observed 
galaxy, with no degree of subjectivity and little sen- 
sitivity to orientation effects. Also, given the TuUy- 
Fisher relation, it is clear that Vd in eqlQ)) can be 
replaced for total magnitude on a given band, for 
cases where the rotation velocity is not available. 

We now study the expected scaling between A 
and other observables of a disk galaxy, for example 
colour, gas fraction, star formation activity, and de- 
gree of flatness of the disk. The main ingredient will 
be the Toomre parameter. 



Q{r) 



K{r)a(r) 
7rG'E(r) ' 



(5) 



This parameter measures the degree of internal 
gravitational stability of the disk, with K{r) account- 
ing for the disruptive shears induced by the differen- 
tial rotation, a{r) modeling the thermal pressure of a 
component, both in competition with S(r)G, the self 
gravity of the disk. Values of Q < 1 are interpreted 
as leading to gravitational instability in a 'cold' disk. 
We replace K{r) in eq© with fl{r) = Vd/r, the an- 
gular velocity of the disk. 

The ratio of disk scale height to disk scale length, 
h/Rd, is a measurable characteristic of galaxies 
which it is easy to show, will scale with A. Assuming 
a thin disk, virial equilibrium in the vertical direction 
yields a relation between h and S, 



h = 



27rGS' 



(6) 



We use this relation for h to replace the gas veloc- 
ity dispersion appearing in equation JSJ for a combi- 
nation of h and the surface density. The dependence 
on S is replaced by a dependence on Md, Mh and 
A to get a new expression for the Toomre's stabil- 
ity criterion, which evaluating radial dependences at 
r = Rd yields. 



A 



(7) 



With F = 1/25, evaluating at Q = 1, the stabil- 
ity threshold suggested by self-regulated star forma- 
tion cycles, the ratio h/Rd is obtained as: 
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Fig. 1. Distribution of values of A through eqQl for the 
complete ASSSG sample, with dots showing the binning. 



Fig. 2. h/Rd ratios for galaxies according to their Hubble 
type, for galaxies in the dGKK sample. 



Rd 390A' ^ ' 

Thinking of star formation as determined by a 
Schmidt law, E = TrGSg , we can use a characteristic 
time for the star formation process in the disk Tsf = 
Sg/E, that gives an idea of the time a galaxy can 
support a star formation rate of constant magnitude. 
From the above relation Tsf becomes 

ttG 2n^GRj , 

at fixed disk mass. Alternatively, one could start 
from an empirical Schmidt law of power n, a scaling 
Tsf oc A" will always result, which does not alter 
our conclusions provided n > 1, a general trait of 
Schmidt laws found in the literature. 

Large values of A will correspond to long star 
formation periods; for this case we expect to have 
galaxies with young populations, looking relatively 
blue and having the high gas fractions of late type 
disks. Galaxies with low A will have short star for- 
mation periods; for them we will hardly see young 
stars and they will look red and gas poor, as is the 
case for early types. A scatter on this trend will 
be introduced by variations in total mass, together 
with a systematic reduction in Tsf in going to larger 
masses (c.f. equation |3, in accordance with larger 
average galactic masses found at earlier types. 

3. DATA AND MODEL COMPARISONS 

As mentioned in the preceding section, we expect 
a correlation between the A parameter and the Hub- 
ble type, since the analysis of section 2 leads one to 
suspect that it is this what primarily determines the 
properties of the galaxy which are taken into account 
for its classification. We now test these expectations 
by direct comparison against two galaxy samples: 



the atlas of Courteau (1996) and the sample from 
Kregel et. al (2005), henceforth ASSSG and dGKK, 
respectively. 

In figure 1 we present the distribution of the 
ASSSG sample as a function of A, estimated through 
eq(0J. For the entire sample, the mean value is 
< A >= 0.0645. In general, the bulk of a galac- 
tic population moves to larger values of A for later 
types, the mean value for each group is: < X >si,= 
0.0574, < A >sbc^ 0.0649, < A >sc= 0.0689. 

We see in figure 1 an approximately Gaussian 
shape, with a considerable extension towards the 
larger values of A. We do not explore the details 
of this distribution further, as the ASSSG sample 
is not intended to be complete or statistically unbi- 
ased in any sense. However, use of eqQ in one such 
sample should prove useful in obtaining an empirical 
distribution of A parameters, to compare for exam- 
ple, against predictions from cosmological models. 

The analysis of the distribution of h/Rd ratios 
was made with the dGKK sample. It contains galax- 
ies with more morphological types; from Sa to Scd 
galaxies. Figure 2 shows early galaxies presenting 
varied values for the h/Rd ratio, although the means 
show a clear diminishing trend in going to latter 
types. Late type galaxies at the Sc end clearly show 
on average much thinner disks than the earlier types, 
as confirmed by several authors over the last few 
years. 

The complementary plot to figure 2 is given in 
figure 3, where instead of presenting h/Rd vs A, we 
use a logarithmic plot to better detect the presence of 
a relation between these two parameters of the type 
the dimensional analysis of section 2 leads us to ex- 
pect, equation®. The line drawn in figure 3 is the 
best fit straight line for the sample, having a slope 
of —0.4 ±0.9, not at odds with value of —1 predicted 
by the simple model. Trends for colour distributions 
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Fig. 3. Logarithmic plot of the relation between A and 
h/Rd, for galaxies in the dGKK sample, the solid line 
shows a linear fit to the data, yielding a slope of —0.4 ± 
0.9. 

and bulge to disk ratios (not shown) arc also qual- 
itatively very similar when plotting against Hubble 
type or A, Hernandez & Cervantes-Sodi (2006). 

We must now check our estimate of A in cases 
where this parameter is in fact known, galaxy forma- 
tion simulations where A acts as an input parameter. 
We use results from two groups, models calculated 
in connection with a study of star formation in disk 
galaxies applied to the Milky Way (Hernandez et al. 
2001), and the models published by van den Bosch 
(2000). 

Both models simulate the cosmological dark mat- 
ter and baryonic mass accretion through a primor- 
dial fluctuation spectrum, and the subsequent, self- 
consistent dynamical and star formation evolution of 
a galaxy. The details of both models vary in numer- 
ical approaches, resolution, time step issues and the 
level of approximation and inclusion of the many dif- 
ferent physical aspects of the complicated problem. 
Figure 4 shows the Xmodei values given as input by 
the above authors to their detailed formation codes, 
with the y-axis showing A from equation 10} applied 
to the final results of the code evolution which fol- 
lowed the model for ^ 13 Gyr. 

It is interesting that the particular features of a 
galactic formation scenario which we are interested 
in treating here. A, final resulting Rd and flat rota- 
tion region Vd, are evidently very well modeled by 
the trivial physics that went into eqQ), across the 
more than 2.5 orders of magnitude covered by the 
masses in the modeled galaxies. This can be under- 
stood by considering that eqQ) is the result of two 
fundamental hypothesis: 
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Fig. 4. A comparison of the input value of A in a series of 
cosmological CDM galactic evolution scenarios, Xmodei, 
and A calculated on the final result of the given sim- 
ulations through eq(|lj. The triangles are models from 
Hernandez et al. (2001) and stars from van den Bosh 
(2000). 

I) an exponentially decreasing surface density profile 
for the disk, and II) a dominant dark halo responsible 
for a flat rotation curve, linked to the disk through 
a Tully-Fisher relation. Whenever the above two 
conditions are met, as is the case for the flnal state 
of all published modeled galaxies and real observed 
systems, basic physics strongly constrains results to 
lie not far from equation 10}. 

To conclude, we have presented an easily appli- 
cable way of estimating the A parameter, a theorist's 
favourite descriptor of a galaxy type, for any real ob- 
served disk galaxy, which is therefore proposed as a 
physically motivated, objective classification param- 
eter for disk galaxies. An extensive development of 
these ideas can be found in Hernandez & Cervantes- 
Sodi (2006). 
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